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ABSTRACT 


The  numerical  methods  relevant  to  the  wake  prediction  of 
submerged  appended  bodies  for  the  purpose  of  predicting  full 
scale  wake  are  surveyed.  The  main  flow  features  are  first 
described  and  the  numerical  techniques  used  to  predict  them 
are  then  discussed.  A  complete  solution  of  the  flows  around 
submerged  appended  bodies  is  beyond  the  reach  of  the  state- 
of-the-art,  both  numerically  and  computerwise.  For  this 
reason,  two  approaches  are  recommended  to  be  pursued 
simultaneously.  One  approach  is  for  inmediate  applications 
and  the  other  approach  for  long-range  research.  For  the 
first  approach,  the  momentum  integral  method  is  recommended 
due  to  its  relative  simplicity  and  its  flexibility  to  adapt 
to  experimental  data.  For  the  second  approach,  a  numerical 
procedure  no  less  sophisticated  than  the  one  by  solving 
the  partially  parabolic  Navier-Stokes  equations  with  the  k-{ 

turbulence  model  should  be  considered.  ,  ,  ? ' 

<  Ka/sa 


ADMINISTRATIVE  INFORMATION 

This  investigation  was  authorized  and  funded  by  the  Office  of  Naval  Technology 
Ships  Submarines  and  Boats  Exploratory  Development  Program,  under  the  management 
of  the  Naval  Sea  Systems  Command,  Code  05R,  Program  Element  6254 3N,  Task  Area 
421-252,  Work  Unit  number  1-1506-202-11. 

INTRODUCTION 

In  order  to  design  a  propeller  for  a  submerged  vehicle,  it  is  essential  to 
know  the  velocity  field  (i.e.  wake)  in  the  propeller  plane.  One  of  the  goals  of 
the  application  of  computational  fluid  dynamics  to  ship  design  is  to  develop 
numerical  techniques  such  that  the  wake  in  the  propeller  plane  can  be  predicted. 
Once  an  analytical  predition  method  has  been  developed,  model  tests  for  the 
purpose  of  obtaining  wake  data  will  be  either  unnecessary  or  would  be  needed 
only  at  the  final  stage  of  design.  Before  such  a  goal  can  be  achieved,  full  scale 
wake  data  must  be  obtained  from  model  wake  data  by  an  appropriate  extrapolation 
method  in  order  to  account  for  Reynolds  number  effect. 


In  the  case  of  submerged  appended  bodies,  numerical  techniques  have  been 
developed  for  the  prediction  of  the  wake  of  bodies  of  revolution,  but  numerical 
techniques  are  not  available  to  predict  the  wake  behind  appended  bodies.  Thus, 
currently  the  full  scale  wake  of  a  submerged  appended  body  is  obtained  by  the 
superposition  of  the  p red i ted  full  scale  wake  of  the  bare  body  and  a  correction 
due  to  the  presence  of  the  appendages.  Since  the  correction  is  based  on  model 
test  data,  the  Reynolds  number  effect  due  to  the  presence  of  the  appendages  has 
been  completely  neglected.  Hie  development  of  numerical  wake  prediction  methods 
which  recognize  the  presence  of  the  appendages  will  undoubtedly  inprove  the 
accuracy  of  the  predictecd  full  scale  wake.  The  objective  of  this  report  is  to 
survey  available  numerical  techniques  which  are  relevant  to  the  wake  prediction  of 
submerged  appended  bodies  and  to  recommend  approaches  which  appear  promising. 

The  wake  predicted  in  the  absence  of  the  propulsor  (i.e.  nominal  wake) 
is  modified  by  the  influence  of  an  operating  propeller.  This  report  will  not 
include  a  survey  of  the  techniques  used  to  account  for  wake/prcpulsor  interactions; 
only  numerical  techniques  which  are  relevant  to  the  prediction  of  the  nominal 
wakes  of  submerged  appended  bodies  will  be  surveyed. 

As  it  will  became  evident  later,  the  flow  around  an  appended  body  has  not 
been  solved  successfully  for  the  application  described  above.  Therefore,  it  is 
not  possible  to  give  a  critical  review  of  the  numerical  techniques  available 
to  solve  the  problem.  Instead,  the  numerical  techniques  which  have  been  used 
in  other  research  fields,  such  as  aerodynamics  and  turbcmachinery,  to  solve  the 
flow  problems  relevant  to  the  one  considered  here  will  be  surveyed  and 
introduced. 
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At  present,  limitation  of  both  the  numerical  technique  and  computer 
capability  make  it  unfeasible  to  attempt  to  compute  the  wake  fields  of  appended 
bodies  by  solving  the  full  Navier-Stokes  equations.  It  is  therefore  necessary  to 
sinplify  the  governing  equations  and,  consequently,  it  is  also  necessary  to  select 
a  turbulence  closure  model  appropriate  to  the  problem  considered.  To  do  both, 
it  is  essential  to  understand  the  physics  of  the  problem  considered.  Thus,  th< 
main  flow  characteristics  of  unappended  bodies  will  be  discussed  first  and  ther 
the  major  complications  of  the  flow  field  due  to  the  presence  of  the  appendage 
will  be  described.  Ihe  discussion  of  various  numerical  methods  and  turbulence 
closure  models  commonly  used  in  turbulent  boundary  layer  calculations  will  follow. 
Advantages  and  disadvantages  of  these  methods  and  models  will  be  mentioned  based 
on  comparisons  of  numerical  methods  that  took  place  in  several  international 
conferences.  This  will  be  followed  by  a  review  of  the  specific  numerical 
techniques  used  in  the  treatment  of  the  main  flow  features  relevant  to  appended 
bodies.  Based  on  the  above  discussion,  recommendation  of  promising  approaches 
will  be  made. 

In  addition  to  the  ninety  two  references  given  for  the  above  discussion,  an 
alphabetical  bibliography  of  the  literature  surveyed,  which  includes  the  ninety 
two  references  mentioned,  is  also  given  in  the  Appendix.  Highlights  of  the 
bibligraphy  are  given  in  Table  A1  of  the  Appendix. 

PHYSICS  OF  THE  PROBLEM 
UNAPPENDED  BODIES  OF  REVOLUTION 

Extensive  measurements  of  mean  velocity,  pressure,  Reynolds  stress  and  skin 
friction  have  been  made  on  various  axisymmetric  bodies  of  revolution  with  emphasis 
on  stern  and  near  wake  regions.  Of  particular  interest  are  the  works  of  the  two 
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research  teams  led  by  Patel1/ 2*  and  Huang3.  The  main  flow  feature  is  that  the 
thickness  of  the  turbulent  boundary  layer  is  in  the  same  order  of  magnitude  as  the 
local  body  radius  in  the  stem  region  at  locations  greater  than,  say,  90  percent 
of  the  total  body  length  (X/L  =  0.9)  for  an  axisymmetric  body,  as  illustrated  in 
Figure  1,  taken  from  Reference  1.  The  main  characteristics  of  a  thick  boundary 
layer  may  be  summarized  as: 

1.  A  large  increase  of  the  pressure  coefficient  fran  the  edge  of  the 
boundary  layer  to  the  wall  as  shown  in  Figure  2.  The  increase  is 
about  70  percent  at  X/L  =  0.95  and  about  100  percent  at  X/L  =  0.99. 

2.  The  mixing  length  in  the  thick  boundary  layer  is  shortened  by  more 
than  50  percent  compared  to  that  in  the  upstream  region  of  the  thin 
boundary  layer  as  shown  in  Figure  3. 

3.  Reynolds  shear  stresses  are  less  intensive  in  the  thick  boundary  layer 
region,  but  spread  out  more.  The  variation  of  the  axial  Reynolds 
shear  stress  is  shewn  in  Figure  4  to  illustrate  this  point. 

Analogous  characteristics  have  also  been  observed  in  the  near  wake  region2. 

Frcm  the  point  of  view  of  seeking  numerical  solutions  in  the  thick  turbulent 
boundary  layer,  characteristic  (1)  implies  that  the  governing  equations  based  on 
the  standard  thin  boundary  layer  approximation  require  modifications  to  account  for 
the  pressure  variation  across  the  bounday  layer,  and  characteristics  (2)  and  (3) 
imply  that  the  turbulence  closure  models,  in  particular  the  eddy-vi scosity  model 
which  is  based  on  the  mixing  length  concept,  must  also  be  modified  accordingly. 

There  are  several  mechanisms  which  cause  the  thick  boundary  layer.  Among 
them  are  the  adverse  pressure  gradient  in  the  streamwise  direction,  the  divergence 
of  mean-flow  streamlines  in  planes  normal  to  the  body  surface  and  the  transverse 
^References  are  listed  on  page  86. 

4 


< 


surface  curvature  effect.  The  first  two  effects  are  the  results  of  a  relatively 
large  streairwise  surface  curvature  in  the  stem  region  of  the  body  of  revolution. 

Hie  effect  of  the  transverse*  surface  curvature  on  the  pressure  variation  across  the 
boundary  layer  is  not  large.  For  example,  as  the  flow  streams  downward  along  the 
surface  of  a  long  slender  cylinder  of  constant  radius  (i.e.  large  transverse 
curvature),  the  boundary  layer  thickness  can  become  several  times  larger  than  the 
radius  of  the  cylinder.  However,  the  rate  of  increase  of  the  displacement  thick¬ 
ness  is  small  and  the  pressure  variation  across  the  boundary  layer  is  not  large 
much  like  the  thin  boundary  layer  of  a  flat  plate.  The  comparison  of  measurements 
and  the  computations  by  Cebeci^  confirms  that  for  this  type  of  flow,  the  thin 
boundary  layer  approximation  is  adequate.  The  above  discussion  leads  to  the  con¬ 
clusion  that  the  main  cause  of  the  thick  boundary  layer  where  the  pressure  varies 
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significantly  across  the  boundary  layer  is  the  convex  longitudinal  surface  curvature 
in  the  stem  region  where  the  streamlines  converge.  As  the  flow  streams  along  the 
curved  surface,  a  large  pressure  gradient  is  generated  in  order  to  balance  the 
centrifugal  force  applied  on  fluid  particles  along  a  curved  streamline. 


APPENDED  BODIES 

Most  of  the  theoretical  and  experimental  investigations  of  flews  around 
appended  bodies  were  carried  out  in  relation  to  the  wing/fuselage  junction  flews 
in  aircrafts  and  the  blade/end-wall  flews  in  turbomachines.  These  two  areas  of 
research  are  the  main  sources  of  this  survey  in  the  attempt  to  understand  the 
physical  phenomena  involved  in  flows  around  appended  bodies  and  to  search  for 
promising  numerical  approaches  to  predict  them. 

The  complexity  of  the  flow  field  around  an  appended  body  is  shewn 
schematically  in  Figure  5  taken  from  Reference  5.  The  resulting  wake  shewn  in 


Transverse  is  used  here  to  mean:  transverse  to  the  streamwise  flew  direction. 
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Figure  6  indicates  the  importance  of  the  interference  effect  between  the  appendage 
and  the  body.  A  further  indication  of  the  inportance  of  the  interference  effect 
between  a  body  of  revolution  and  a  tail  fin  has  been  demonstrated  in  the  wake 
measurements  by  Sevik^ 

The  flew  field  of  the  appendage/body  junction  as  shown  in  Figure  5  can 
be  divided  into  three  regions: 

1.  The  longitudinal  flcw-reversal  region:  This  is  the  region  between  the 
singular  separation  point  and  the  leading  edge  of  the  appendage  (strut 
stagnation  line).  The  horseshoe  vortex  is  generated  in  this  region. 

2.  The  transverse  flcw-reversal  region:  This  region  is  the  passage  of  the 
horseshoe  vortex  along  the  two  sides  of  the  appendage/body  junction  and 
is  downstream  of  the  longitudinal  flcw-reversal  region. 

3.  The  crossflcw  region:  This  is  the  region  outside  of  the  above  two 
regions  that  are  enclosed  by  the  line  of  ordinary  separation,  or  what 
is  also  called  the  primary  separation  line. 

These  flow  features  have  been  observed  by  flow  visualization  techniques^ f 7-10 
and  can  be  explained  as  follows. 

Suppose  the  appendage  (a  wing  or  a  cylinder)  is  mounted  on  a  flat  plate 
or  on  an  axisymmetric  body  of  revolution,  then  the  boundary  layer  far  upstream  of 
the  appendage  is  two-dimensional.  The  two-dimensionality  of  the  flow  will  be 
destroyed  by  the  pressure  force  generated  by  the  presence  of  the  appendage.  This 
pressure  force  tends  to  push  the  fluid  away  from  the  appendage .  Since  the  fluid 
particles  near  the  wall  of  the  body  have  been  slowed  down  by  the  wall  friction, 
the  inertial  force  near  the  wall  is  smaller  than  the  inertial  force  away  from  the 
wall  and  is  more  easily  overcome  by  the  pressure  force.  Thus,  the  interaction 
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between  the  inertial  and  the  pressure  forces  causes  the  flow  to  reverse  in  the 
longitudinal  direction  in  region  (1)  as  shown  in  Figure  5,  generating  the  horseshoe 
vortex.  The  vorticity  line  of  the  horseshoe  vortex  then  bends  around  the  appendage 
like  a  necklace,  resulting  in  streanwise  vortices  in  region  (2)  with  their 
directions  as  shewn  in  Figure  5.  There  are  two  types  of  crossflows  in  region  (3), 
one  is  the  non-crossover  type  as  shewn  in  Figure  7  and  the  other  is  the  crossover 
type  as  shown  in  Figure  8.  The  non-crossover  crossflow  occurs  in  the  upstream  and 
the  downstream  position  away  from  the  inflection  point  of  the  S-shaped  streamline 
as  shewn  in  Figure  5  while  the  crossover  crossflow  occurs  in  the  neighborhood  of 
the  inflection  point  where  the  direction  of  the  curvature  of  the  free-stream 
streamline  is  reversed.  The  direction  of  the  crossflow  upstream  of  the  inflection 
point  is  the  same  as  that  shewn  in  Figure  7,  but  the  direction  is  reversed 
downstream  of  the  inflection  point  since  the  direction  of  the  curvature  of  the 
freestream  streamline  is  reversed  resulting  in  the  reversal  of  the  direction  of 
the  centrifugal  force.  Both  types  of  crossflows  are  pressure-driven.  Consequently, 
rapid  changes  in  the  crossflow  profile  take  place  near  the  wall  and  little  changes 
occur  near  the  edge  of  the  boundary  layer  where  the  pressure  force  is  approximately 
balanced  by  the  centrifugal  force  generated  by  tie  curvature  of  the  streamline. 

For  a  further  discussion  of  both  types  of  crossflows,  one  is  referred  to  the  papers 
by  Klinksiek  and  Piercell  and  by  Prahladl2. 

The  laminar  horseshoe  vortex  has  been  observed  to  form  in  systems  of  2, 

4,  or  6  vortices  as  shown  schematically  in  Figure  9  (taken  from  Reference  7). 
However,  two  vortices  are  predominant  in  the  turbulent  boundary  layer  and  further¬ 
more,  these  two  vortices  tend  to  merge  into  one  as  shewn  by  Peake  et  al®  in  Figure 
10.  This  is  further  confirmed  by  the  measurement  of  the  velociy  vector  dewastream 
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of  the  leading  edge  of  the  wing  by  Shabaka  and  Bradshaw13  shown  in  Figures  11a 
and  lib.  The  implication  of  this  observation  is  that,  if  the  detailed  flow 
features  in  region  (1)  are  not  of  major  concern,  the  flow  field  can  be  assumed 
as  being  generated  by  a  single  horseshoe  vortex.  This  assumption  is  important 
in  modelling  the  streamwise  vortex  motion  in  region  (2)  when  momentum  integral 
methods  are  used. 

The  separation  point  at  the  leading  edge  of  region  (1)  is  a  singular  point 
with  zero  skin  friction1*).  But  the  skin  friction  is  not  zero  along  the  primary 
separation  line1*);  rather  it  is  a  line  along  which  the  wall  limiting  streamlines 
from  region  (3)  and  those  from  regions  (1)  and  (2)  converge.  Further  discussion 
of  the  type  and  of  the  nature  of  flow  separations  can  be  found  in  References  14-17. 

Another  feature  of  turbulent  boundary  layers  with  flow  separations  or  with 
large  crossflows  is  that  the  direction  of  the  velocity  gradient  vector  differs 
significantly  from  that  of  the  shear  stress  vector  even  near  the  wall.  This 
anisotropic  property  of  the  shear  stress  has  been  observed  by  several  investigators 
under  different  experimental  conditions1*)-2*).  The  inplication  of  this  property  is 
that  the  sinple  eddy-viscosity  and  the  mixing-length  models,  which  automatically 
lead  to  the  coincidence  of  these  two  vectors,  may  not  be  adequate  turbulence 
closure  models  for  this  type  of  problem. 

The  size  of  the  horseshoe  vortex  region  depends  on  the  geometry  of  the 
appendage*) , 21 ,  and  more  importantly  on  the  thickness  of  the  incoming  turbulent 
boundary  layer  near  the  location  of  the  singular  separation  point*)*21.  Thus, 
referring  to  a  body  of  revolution  (see  Figure  1),  if  a  tail  fin  is  placed  at  95 
percent  length  of  the  body  where  the  thickness  of  the  boundary  layer  is  approxi¬ 
mately  equal  to  the  local  radius  of  the  body,  then  the  size  of  the  horseshoe  vortex 
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region  is  approximately  in  the  order  of  the  local  radius  of  the  body.  On  the 
other  hand,  if  a  fin  is  placed  close  to  the  nose  of  the  body  then,  the  horseshoe 
vortex  region  generated  will  be  rather  small  since  the  thickness  of  the  bounday 
layer  at  that  location  is  very  thin. 

NUMERICAL  METHODS 

Several  papers  have  discussed  the  state-of-the-art  and  the  outlook  of 
computational  fluid  dynamics.  In  the  1980  symposium  on  the  "Computation  of  Viscous  - 
Inviscid  Interaction",  Le  Balleur22  divided  the  numerical  approximations  used  in  the 
ccrputation  of  laminar  and  turbulent  flows  into  four  levels  and  indicated  that  the 
state-of-the-art  is  no  higher  than  the  middle  of  the  third  level  as  indicated  by 
symbol  x  in  Table  1.  At  this  state-of-the-art,  only  the  problems  which  are  suitable 
for  a  strong  viscous/inviscid  interaction  based  on  the  time-averaged  Navier-Stokes 
equations  with  the  thin  boundary  layer  approximation  can  be  solved  with  confidence. 
Kline,  Ferziger  and  Johnston23  and  Chapman24  discussed  computational  fluid  dynamics 
(1978  and  1979  respectively) .  Summaries  of  these  outlooks  are  presented  in  Tables 
2  and  3. 

To  give  seme  indications  as  to  when  the  problem  of  the  flows  around  appended 
bodies  can  be  solved  using  the  full  Navier-Stokes  equations,  a  projection  by 
Chapman24  on  a  similar  problem  in  aerodynamics  is  worthwhile  mentioning.  According 
to  him  the  practical  3D  wing/body  problem  probably  can  be  solved  using  the  full 
*  Navier-Stokes  equations  by  1990’s,  but  still  at  a  prohibitive  cost  and,  ir.  addition, 

'  a  more  than  two  orders  of  magnitude  larger  computer  than  what  is  available  today 

would  be  required. 

The  above  discussion  indicates  that  at  present  the  attempt  to  predict  the  wake 
of  appended  bodies  by  solving  the  full  Navier-Stokes  equations  is  unrealistic. 
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Therefore,  the  survey  of  the  numerical  methods  adequate  for  wake  prediction  will 
be  confined  to  those  associated  with  solving  the  time-averaged  Navier-Stokes 
equations  with  various  degrees  of  sinplifications. 

Because  of  the  above  restriction,  several  interesting  numerical  methods  will 
not  be  discussed.  Among  than  are  the  spectral  method25,  the  large  eddy  simula¬ 
tion  with  a  subgrid  closure26  and  the  discrete  vortex  method27,28.  of  particular 
interest  is  the  recent  development  of  the  discrete  vortex  method  which  appears  to 
be  rather  attractive  to  simulate  separated  flows  and  wakes.  For  a  quite  different 
reason,  the  finite  elenent  method  will  also  not  be  discussed.  It  is  the  author's 
opinion  that  the  choice  between  a  finite  element  method  and  a  finite  differential 
method  is  more  of  a  matter  of  personal  preference;  furthermore,  the  finite  element 
method  is  used  infrequently  in  the  computation  of  turbulent  boundary  layers.  Some 
fundamental  aspects  of  the  finite  elenent  method  in  fluid  mechanics  are  discussed 
by  Shen29. 

Thus,  the  numerical  approaches  to  be  surveyed  are  finite  differential  methods 
and  momentum  integral  methods  based  on  the  simplified  Navier-Stokes  equations.  When 
the  thin  boundary  layer  approximation  is  used,  or  even  when  this  approximation  is 
relaxed,  the  pressure  distribution  needs  to  be  provided  in  order  to  initiate  boundary 
layer  calculations.  This  pressure  distribution  is  normally  provided  by  the  potential 
flow  calculation.  For  this  reason,  calculation  methods  for  the  potential  flew 
will  also  be  discussed. 

FINITE  DIFFERENTIAL  METHODS 

Depending  on  approximations  made,  the  governing  equations  can  be  divided  into 
four  types  as  shown  in  Table  4.  The  mathematical  classification  of  the  typ>es  of 
the  governing  equations  can  be  found  in  Wang's  papier 34. 


Fran  a  computational  point  of  view,  the  distinctions  among  these  four  types  of 
governing  equations  are  important.  For  the  hyperbolic  and  the  parabolic  type 
of  equations,  the  pressure  field  is  assumed  known  and  the  velocity  components  can 
be  computed  by  marching  downstream.  No  iteration  of  solutions  is  needed  in  this 
calculation  procedure.  For  the  elliptic  type  of  equations,  both  the  pressure  and 
the  velocity  components  are  unknowns  and  must  be  obtained  by  iteration  of  solutions, 
thus  resulting  in  large  demands  on  the  computer  memory  and  the  computer  cost. 

Itie  complexity  of  the  calculation  procedure  for  the  partially  parabolic  type  is 
intermediate  between  the  above  two  extremes.  The  pressure  field  of  the  entire 
computational  domain  is  assumed  known  and  the  velocity  components  are  then  computed 
by  marching  downstream.  Once  these  are  done,  a  new  pressure  field  is  calculated 
from  than  and  the  calculational  procedure  is  repeated  until  a  convergence  of  the 
solution  is  obtained.  Thus,  unlike  the  elliptic  type  of  equations,  only  a  three- 
dimensional  memory  of  the  pressure  is  required  by  the  partially  parabolic  type. 

The  calculation  procedure  for  the  partially  parabolic  type  was  first  applied 
by  Pratap  and  Spalding^2  to  the  flow  in  curved  ducts.  The  results  agree  with  the 
experimental  results  better  than  those  obtained  by  the  parabolic  type  of  equations. 

A  further  comparison  of  these  two  types  of  equations  has  been  discussed  by  Spalding 
and  his  ccworkers35.  Because  the  pressure  field  of  the  entire  computational  domain, 
including  the  normal  pressure  gradient,  is  considered,  the  use  of  the  calculation 
procedure  of  the  partially  parabolic  type  of  equations  is  worthy  of  a  careful 
consideration  for  the  prediction  of  the  flows  around  appended  bodies. 

To  close  the  time-averaged  Navier-Stokes  equations  such  that  the  number  of 
equations  and  the  number  of  unknowns  are  equal,  a  turbulence  closure  model  is  needed 
to  connect  the  Reynolds  shear  stress  and  the  mean  field.  Turbulence  closure  models 
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being  used  at  present  can  be  divided,  according  to  the  order  of  increasing 
complexity,  into:  zero-equation  model,  one-equation  model,  two-equation  model, 
stress-equation  model  and  large-eddy  simulation  with  a  subgrid  scale  model. 
Detailed  discussions  of  these  models  can  be  found  in  the  review  papers  by 
Reynolds^  and  Reynolds  and  Cebeci37.  The  experience  with  numerical  computations 
indicates  that  the  use  of  the  stress-equation  model  where  each  component  of 
the  six  Reynolds  shear  stresses  is  solved  by  a  partial  differential  equation 
is  uneconomical  at  present  in  the  sense  that  the  improvement  in  the  quality  of 
solution  is  insignificant  despite  a  considerable  increase  in  the  amount  of  the 
computational  effort.  The  use  of  the  large-eddy-simulation  model  is  even  more 
complicated.  Thus,  the  use  of  the  zero  equation,  one-equation  and  two-equation 
model  appears  to  be  worthy  of  further  explorations  for  the  purpose  considered 
here,  despite  the  deficiencies  of  these  models. 

The  two  most  popular  versions  of  the  zero-equation  model  are  the  eddy- 
viscosity  model  and  the  mixing-length  model.  The  basic  assumption  for  both  is 
that  the  local  production  of  the  turbulence  energy  is  balanced  by  the  dissipation 
of  the  turbulence  energy.  Consequently,  the  Reynolds  shear  stress  is  related 
to  the  mean  velocity  field  by  an  algebraic  equation.  This  model  has  been 
fully  exploited  by  Cebeci37.  por  complex  turbulent  flews,  such  as  flews  around 
appended  bodies,  the  assumption  that  the  local  turbulence  energy  production  is 
equal  to  its  dissipation  is  not  valid,  and  the  advection  and  diffusion  of 
the  turbulence  energy  must  be  included  in  the  energy  budget.  Therefore,  a  higher- 
order  turbulence  closure  model  may  be  needed.  The  most  developed  one-equation 
model  is  that  proposed  by  Bradshaw  et  al^O  where  the  equation  for  the  turbulent 
kinetic  energy  is  solved  for  the  connection  between  the  Reynolds  shear  stress 


and  the  mean  field.  This  model  has  an  excellent  performance  in  the  1968  Stanford 
conference^ .  For  a  two-equation  model,  not  only  the  equation  for  the  turbulent 
kinetic  energy  but  also  the  equation  for  the  dissipation,  or  a  characteritsic 
turbulent  length  scale  is  solved.  The  construction  of  two-equation  models  is 
currently  an  active  research  area.  A  critical  examination  of  four  popular  two- 
equation  models  has  been  given  in  Reference  39. 

As  canton  in  the  numerical  calculation  of  turbulent  boundary  layers,  the 
choice  of  which  turbulence  closure  model  is  most  appropriate  for  the  flow  around 
appended  bodies  is  not  trivial  nor  clear-cut  even  when  the  physics  of  the  problem 
is  relatively  well  understood.  In  a  recent  calculation^  of  3D  boundary  layers 
involving  flaw  separations  on  practical  wing  configurations,  despite  the  fact  that 
the  sinple  eddy-viscosity  model  is  inadequate  for  the  physics  of  the  problem,  the 
results  obtained  by  more  sophisticated  models  such  as  one-equation  and  two- 
equation  models  dc  not  seem  to  significantly  improve  those  obtained  by  using  a 
sinple  eddy-viscosity  model.  The  conclusion  is  that  a  simple  turbulence  model, 
such  as  the  eddy-viscosity  model  or  Bradshaw's  one-equation  model 30  should 
be  preferred  subject  to  modifications  according  to  the  specific  physics  of  the 
problem  under  consideration. 

MOMENTUM  INTEGRAL  METHODS 

One  of  the  major  difficulties  in  the  calculation  of  turbulent  boundary  layers 
is  the  specification  of  the  Reynolds  shear  stress  everywhere  in  the  boundary  layer. 
One  of  the  main  motivations  of  momentum  integral  methods  is  to  avoid  the  need  to 


specify  the  Reynolds  shear  stress  everywhere  in  the  boundary  layer  except  on  the 
wall.  This  is  done  by  integrating  the  governing  equations  from  the  wall  to  the 
edge  of  the  boundary  layer  such  that  only  the  skin  friction  on  the  wall  needs  to 
be  specified,  which  is  then  obtained  from  empirical  data. 

In  addition  to  the  eirpirical  skin  friction  formula  and  the  mean  velocity 
profile  including  both  the  streamwise  and  the  transverse  (i.e.  crossflow)  direc¬ 
tion,  an  auxiliary  equation  is  needed.  Depending  on  which  of  the  following  three 
equations  -  the  entrainment  equation4!,  the  energy  equation42  and  the  moment  of 
momentum  equation^  is  chosen  as  the  auxiliary  equation,  momentum  integral  methods 
are  divided  into  three  families44.  since  the  entrainment  equation  appear  to  be 
most  thoroughly  developed,  the  use  of  it  as  the  auxiliary  equation  is  recommended. 

In  the  recent  past,  momentum  integral  methods  have  received  severe  criticisms 
from  several  people,  including  Wheeler  and  Johnston45,  Landweber4^  and  Landweber 
and  Patel4?.  Hie  main  criticism  is  in  the  use  of  mean  velocity  profiles.  Of 
particular  concern  is  the  fact  that  the  crossflow  velocity  profile  varies  so  much 
according  to  different  flow  problems  considered  that  it  is  impossible  to  find  a 
universal  crossflow  velocity  profile  which  is  applicable  to  a  wide  class  of 
problems.  However,  if  the  interest  is  confined  to  a  specific  problem  where  drastic 
changes  in  flow  characteristics  are  not  expected,  then  the  ground  for  the  above 
criticism  disappears.  The  flows  around  appended  bodies  considered  here  belong  to 
this  class  of  problems  and  momentum  integral  methods  deserve  a  serious  consideration 
because  of  their  simplicity  of  computation  and  their  flexibility  to  adjust  to 
experimental  data. 

For  momentum  integral  methods,  the  streamline  coordinate  system  appears  to 
be  most  adequate.  As  mentioned  above,  empirical  formulas  for  the  skin  friction, 
the  mean  velocity  profiles  and  the  entrainment  function  are  needed  in  these  methods. 
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These  formulas  are  mostly  obtained  from  the  measurements  of  various  two-dimensional 
turbulent  boundary  layers.  If  the  streamline  coordinate  system  is  used,  the 
turbulent  boundary  layer  in  the  streairwise  direction  of  the  free  stream  can  be 
regarded  as  "two-dimensional"  and  those  formulas  mentioned  above  can  be  applied 
directly  in  the  streanwise  direction.  Then  the  skin  friction  and  the  entrainment 
function  in  the  crosswise  direction  can  be  determined  once  the  limiting  crossflow 
angle  (i.e.  the  angle  between  the  wall  limiting  viscous  flow  direction  and  the 
free  stream  streanwise  direction)  is  computed  from  a  given  crossflow  velocity 
profile.  This  approximation  essentially  reduces  the  additional  information 
required  by  three-dimensional  turbulent  boundary  layers  to  the  specification  of 
a  crossflow  velocity  profile.  The  justification  of  the  use  of  a  two-dimensional 
turbulent  boundary  layer  mean  velocity  profile  family  to  represent  a  three- 
dimensional  turbulent  boundary  layer  has  been  critically  discussed  and  confirmed 
by  Cuirpsty4**. 

POTENTIAL  FLOW  METHODS 

The  progress  in  the  calculation  of  the  potential  flow  around  an  arbitrary 
configuration  such  as  an  aircraft  has  been  rapid  and  successful  in  recent  years. 

The  main  driving  forces  of  the  success  are  the  use  of  the  panel  method  to  describe 
the  surface  of  a  complicated  geometry  and  the  rapid  increase  in  recent  years  of 
the  computer  capacity  which  makes  the  use  of  the  panel  method  feasible.  For  the 
state-of-the-art  of  potential  flow  methods,  one  is  referred  to  the  papers  by  Hess49 
and  by  Carmichael  and  Erickson50, 

The  classification  of  potential  flow  methods  is  non-trivial  and  often 
mathematically  inconsistent.  It  suffices  for  the  purpose  here  to  divide  them  into 
zero-order  methods  and  higher-order  methods.  In  a  zero-order  method,  the  panel 
element  is  flat  and  the  source  distribution  over  it  is  a  constant  while  in  a 


higher-order  method,  either  the  panel  element  is  not  flat  or  the  source  distribution 
is  not  a  constant,  or  both.  The  XYZPF  computer  program  used  at  DINSRDC  is  a  zero- 
order  method  based  on  a  theory  developed  by  Hess  and  Smith  (cf.  Hess4^).  In 
this  method,  the  panel  element  is  a  flat  quadrilateral  and  the  source  strength  is 
constant.  The  theory  developed  by  Webster^l  is  a  higher-order  method  where  the 
source  strength  over  it  is  distributed  as  a  linear  function  of  position  over  the 
panel  element  which  is  a  flat  triangle.  The  most  developed  higher-order  method 
is  the  PAN  AIR  computer  program  discussed  in  Reference  50. 

The  main  short  comings  of  the  zero-order  method  such  as  the  one  used  in  the 
XYZPF  program,  are  that  the  source  distribution  is  discontinuous  between  panels 
resulting  spurious  vortex  sheets  and  that  leakage  occurs  between  panels  because 
flat  trapezoids  cannot  properly  cover  a  strongly  curved  surface,  such  as  the  nose 
or  the  tail  portion  of  an  ellipsoid.  It  is  possible,  that  by  increasing  the  number 
of  panels  the  difficulties  mentioned  above  may  disappear,  at  least  partially. 
However,  for  a  configuation  with  a  complicated  geometry  such  as  the  one  under 
consideration,  the  number  of  panels  may  not  be  allowed  to  increase  to  a  desirable 
level  because  of  the  limitation  in  the  computer  storage.  Thus,  for  the  problem 
considered  here,  a  higher-order  method  is  reccrmended . 

There  is  a  further  consideration  which  also  favors  the  use  of  a  higher- 
order  method.  When  a  simplified  version  of  the  Navier-Stokes  equations,  such  as 
the  parabolic  typ>e,  is  used  to  solve  coup  lex  turbulent  flows,  an  adequate 
scheme  of  the  viscous/inviscid  interaction  is  essential  for  obtaining  adequate 
'  solutions.  As  discussed  above,  the  complex  flow  field  around  an  appended  body 

includes  flew  reversals  in  the  streamwise  and  the  transverse  direction,  both  driven 
by  pressure  forces.  Since  the  ellipticity  of  tine  flow  has  been  eliminated  by  the 
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parabolization  of  the  governing  equations,  the  only  mechanism  by  which  downstream 
conditions  can  influence  upstream  stations  is  through  the  pressure  distribution 
which  is  calculated  by  the  potential  flew  method.  Thus,  accurate  potential  flow 
calculations  are  crucial  in  obtaining  viscous  flow  solutions  and  a  higher-order 
calculation  method  for  the  potential  flow  appears  to  be  more  adequate. 

COMPARISON  OF  NUMERICAL  METHODS 

In  several  international  conferences,  some  test  cases  with  experimental  data 
were  set  up  and  researchers  with  numerical  codes  were  invited  to  perform 
confutations  on  the  test  cases;  the  computational  results  were  then  compared 
and  the  numerical  methods  evaluated.  Table  5  provides  a  list  of  six  major 
conferences  with  brief  comments  on  the  findings;  the  emphasis  here  is  on  the 
comparison  between  finite  differential  methods  and  momentum  integral  methods. 

Of  particular  interest  is  the  confutation  of  the  flew  induced  on  a  flat 
plate  by  a  circular  cylinder  based  on  the  experimental  data  of  Bast  and  Hoxeyl°. 

It  is  considered  a  reasonable  goal  of  the  current  state-of-the-art  of  three- 
dimensional  turbulent  boundary  layer  calculation  methods  to  predict  the  separation 
line.  Yet,  none  of  the  methods  presented  by  the  participants  was  capable  of  doing 
this  in  1975  (cf.  East52),  itie  calculation  of  one  of  the  test  cases  in  the 
1980-81  Stanford  Conference^  on  the  junction  flew  in  an  idealized  wing/body 
combination  based  on  the  experiment  by  Shabaka  and  Bradshawl3  (cf.  Figures  11a 
and  lib)  is  also  of  interest.  However,  the  result  of  evaluation  is  not  available 
at  the  time  of  this  writing. 

A  general  conclusion  on  the  comparison  of  finite  differential  methods  and 
momentum  integral  methods  can  be  derived  from  the  results  of  the  conferences  listed 
in  Table  5.  Despite  the  potential  advantage  of  finite  differential  methods  in  the 


coiputation  of  more  ccnplex  flew  problems,  such  an  advantage  has  not  been 
demonstrated  in  general  validity  at  present;  momentum  integral  methods  continue 
to  perform  as  well  as  finite  differential  methods  in  all  classes  of  flow  problems 
tested  including  those  with  flew  separations.  As  concluded  in  the  preliminary 
Evaluation  Conti ttee  Report  of  the  1980-81  Stanford  Conference  based  on  the 
various  test  problems  considered,  momentum  integral  methods  continue  to  perform 
adequately  and  for  engineering  purposes  are  sufficient  and  sometimes  preferable. 

NUMERICAL  TREATMENT  OF  MAIN  FLOW  FEATURES  RELEVANT  TO  APPENDED  BODIES 

If  tail  fins  are  attached  to  the  axisymmetric  body  of  revolution  shewn  in 
Figure  1  at  a  location  approximately  95  percent  of  the  total  length  downstream 
of  the  nose  of  the  body,  then  the  tail  fins  will  be  irmersed  in  the  thick 
turbulent  boundary  layer  as  mentioned  earlier.  The  main  flow  features  in  the 
boundary  layer  and  the  wake  of  the  tail-fin/body  combination  can  be  summarized 
as:  thick  turbulent  boundary  layer,  two-dimensional  and  three-dimensional  flow 
separation,  secondary  flew  and  wake.  The  wakes  of  tail  fins  will  interact 
with  the  bounday  layer  of  the  body.  Hcwever,  this  interaction  is  unlikely  to  be 
important  since  the  bounday  layers  on  tail  fins  are  expected  to  be  small  compared 
to  that  on  the  body. 

THICK  TURBULENT  BOUNDARY  LAYER 

The  streanwise  vortex  motion  in  the  junction  flew  of  a  fin/body  combination 
depends  on  the  strength  and  the  size  of  the  horseshoe  vortex  generated  in  the 
leading  edge  of  the  fin/body  combination.  Since  the  strength  and  the  size  of 
the  horseshoe  vortex  depends  on  the  thickness  and  the  mean  velocity  profile  of 
the  incident  turbulent  boundary  layer,  an  accurate  prediction  of  the  thick 
turbulent  boundary  layer  is  essential  for  wake  prediction. 
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The  main  difficulty  in  the  calculation  of  the  thick  turbulent  boundary  layer 
arises  due  to  the  approximation  made  that  the  normal  pressure  variation  across  the 
boundry  layer  is  zero  in  order  to  sinplify  the  Navier-Stokes  equations.  Thus,  if 
no  approximation  is  made  the  difficult  does  not  arise,  but  then  the  pressure 
must  be  calculated  simultaneously  with  velocity  carponents.  On  the  other  hand, 
if  the  pressure  is  assumed  to  be  constant  across  the  boundary  layer,  the  equations 
are  simplified  to  the  parabolic  or  the  hyperbolic  type,  but  then,  an  adequate 
viscous/inviscid  interaction  scheme  is  necessary  to  account  for  the  correction 
due  to  the  normal  pressure  variation  across  the  boundary  layer. 

Major  methods  used  for  the  numerical  treatment  of  the  thick  bounday  layer 
are  listed  in  Table  6.  The  order  of  listing  is  according  to  the  decreasing 
degree  of  demand  on  computing  efforts.  Representative  references  are 
also  given  for  each  method. 

As  mentioned  earlier  in  the  experimental  aspect  of  the  thick  turbulent 
boundary  layer,  the  structure  of  turbulence  is  different  from  that  of  a 
flat  plate.  Since  most  of  turbulence  closure  models  are  constructed  based  on  the 
data  obtained  in  flat  plates,  modifications  are  needed  when  applied  to  thick 
turublent  boundary  layers.  In  the  simple  mixing-length  model  used  by  Wang  and 
Huang59,  such  modifications  have  been  made.  The  question  of  whether  modifications 
are  needed  in  two-equation  turbulence  models  does  not  appear  to  have  received 
attention. 

The  approach  adopted  by  Wang  and  Huang  in  making  correction  to  the  mixing 
length  in  the  region  of  the  thick  bounday  layer  is  solely  based  on  experimental 
data.  A  more  rational  approach  has  been  suggested  by  Bradshaw62,63.  as  mentioned 
earlier,  the  thick  boundary  layer  normally  occurs  in  the  region  where  the  surface 
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curvature  along  which  the  flow  travels  is  relatively  large.  In  the  case  of  the 
axisymmetric  body  of  revolution,  this  region  occurs  near  the  stem  where  the 
longitudinal  surface  curvature  is  large.  Bradshaw  preposed  that  the  change  in 
the  structure  of  the  turbulent  boundary  layer  be  accounted  for  by  making  correction 
to  the  dissipation  length  scale  due  to  the  extra  strain  introduced  by  the  effect 
of  the  streamline  curvature.  This  idea  has  been  applied  by  Patel  and  Lee64  to 
the  body  of  revolution  using  a  finite  differential  method  with  Bradshaw's  one- 
equation  model30.  a  similar  application,  but,  using  the  momentum  integral  method 
has  been  done  by  Green,  Weeks  and  Brooman^ . 

A  new  development  to  account  for  the  normal  pressure  variation  across  the 
boundary  layer  using  momentum  integral  methods  is  proposed  by  Le  Balleur61. 

The  departure  of  this  approach  from  the  classical  one  is  that  the  viscous  and 
inviscid  equations  are  solved  simultaneously  in  the  boundary  layer  with  the  boundary 
condition  of  the  inviscid  flow  on  the  wall  determined  by  the  surface  transpiration 

method.  The  ingenuity  of  this  approach  is  that  the  pressure  in  the  boundary 

layer  is  set  equal  to  that  of  the  inviscid  flow  in  the  first-order  approximation 
when  curvature  effects  are  neglected.  Since  at  high  Reynolds  numbers,  the 
experimental  evidence  indicates  that  the  actual  pressure  distribution  is  very 
close  to  that  of  the  inviscid  flow,  this  approach  allows  for  a  variation  in  the 
normal  pressure  across  the  boundary  layer  while  solving  the  same  set  of  the  govern¬ 
ing  equations  used  in  the  thin  boundary  layer  approximation.  The  disadvantage  of 

this  approach  is  that  the  various  boundary  layer  thicknesses  are  defined  based  on 

the  inviscid  flew  velocity  on  the  wall  rather  than  on  the  edge  of  the  boundary 
layer.  Therefore  the  application  of  the  empirical  formulas  which  have  various 
boundary  layer  thicknesses  as  parameters  requires  modifications.  Fhst66  has 
extended  the  theory  to  include  the  second-order  effects  due  to  surface  curvatures 
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in  two-dimensional  problems.  The  extension  to  three-dimensional  problems  has  not 
been  made,  but  it  is  feasible. 

FLOW  SEPARATION 

TVo  types  of  flow  separation  occur  in  the  flow  field  as  shown  in  Figure  5. 

The  separation  point  upstream  of  the  leading  edge  of  the  appendage  is  a  two- 
dimensional  singularity  with  zero  skin  friction;  but  the  separation  line  which 
encircles  the  appendage  like  a  necklace  is  not  a  line  of  singularities,  rather,  it 
is  a  line  along  which  the  limiting  wall  streamlines  converge.  Thus,  numerically, 
the  first  type  of  separation  creates  much  more  difficulty  than  the  second  type. 

Since  the  skin  friction  approaches  zero,  a  natural  criterion  for  the  prediction 
of  the  singular  separation  point  is  to  locate  the  point  where  the  skin  friction  is 
zero.  However,  it  has  been  found  that  numerically  it  is  easier  to  locate  this 
point  approximately  as  where  the  displacement  thickness  increases  sharply.  To 
numerically  march  across  the  singular  separation  line,  the  so  called  inverse  mode 
is  used.  In  the  inverse  mode,  the  displacement  thickness  is  specified  and  the 
pressure  gradient  predicted;  this  is  contrary  to  the  usual  direct  mode  where  the 
pressure  gradient  is  specified  and  the  displacement  thickness  predicted.  A  clear 
illustration  that  the  singularity  in  the  governing  equations  is  removed  by  using 
the  inverse  mode  is  given  by  Cousteix  and  Houdeville67.  Other  applications  of 
the  inverse  mode  to  calculate  separated  flows  using  momentum  integral  methods  are 
discussed  by  Cousteix  et  al68.  East  et  al69,  Whitfield  et  al?0  and  Stock43. 

Similar  applications  using  finite  differential  methods  are  discussed  by  Catherall 
and  Mangier71,  Carter72  and  Formery  and  Delery7^. 

The  paper  by  Formery  and  Delery  is  of  particular  interest.  They  assumed  that 
the  pressure  is  constant  in  the  direction  normal  to  the  wall  but  is  allowed  to 
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vary  in  the  transverse  direction  and  succeeded  in  calculating  the  flow  field  in 
the  vicinity  of  the  separation  line  in  the  experiment  of  East  and  HoxeylO  as  shown 
in  Figure  12.  They  did  not  attempt  to  calculate  the  entire  region  of  the  horseshoe 
vortex  motion  (i.e.  region  (1)  in  the  previous  discussion)  and  probably  would  have 
failed  since  the  elliptic  type  of  the  governing  equations  must  be  used  in  this 
region.  Same  attenpts^^/75  to  calculate  the  flow  field  in  this  region  have 
been  limited  to  low  Reynolds  numbers  in  the  range  of  100-500. 

Since  the  pressure  force  is  the  dominant  driving  force  in  this  region,  the 
governing  equations  chosen  must  allow  for  pressure  variations  in  all  three 
dimensions  or  a  strong  viscous/inviscid  interation  scheme  must  be  used  to  account 
for  pressure  variations  via  potential  solutions  which  are  elliptic.  On  the  other 
hand,  because  the  pressure  is  the  dominant  force,  Reynolds  shear  stresses  are  less 
inportant  and  a  simpler  turbulence  closure  model  may  be  adequate  as  is  suggested 
by  the  calculation  of  Forme ry  and  Delery. 

SECONDARY  FLOW 

Hie  term  secondary  flew  is  used  here  to  indicate  the  flew  field  with  a 

significant  component  of  crossflow  due  to  the  geometry.  This  flow  field  includes 

region  (2)  of  the  junction  flew  with  a  strong  streanwise  vorticity  and  region  (3) 

of  the  large  crossflow  region.  The  latter  region  includes  the  crossover  (or 

S-shaped)  crossflow  in  the  vicinity  where  the  free-stream  streamline  has  a  point 

of  inflection  and  the  non-crossover  crossflow  in  the  upstream  and  the  downstream 

portion  of  the  point  of  inflection.  Unlike  region  (1)  where  the  horseshoe  vortex 

is  generated,  the  flow  reversal  in  the  streanwise  direction  does  not  take  place 
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in  regions  (2)  and  (3).  Calculations  made  by  Pratap  and  Spalding  indicate 
that  a  partially  parabolic  type  of  equations  will  be  adequate  in  these  flew  regions. 


The  experiment  by  Shabaka  and  Bradshawl3  indicates  that  the  direction  of  the 
velocity  gradient  vector  differs  significantly  from  that  of  the  shear  stress  in  the 
junction  flow,  This  implies  that  the  sinple  eddy-viscosity  and  the  mixing-length 
turbulence  models  may  not  be  adequate  in  this  region,  although  the  contrary  is 
indicated  earlier 40.  Numerical  calculations  of  the  junction  flow  in  a  wing/body 
combination  are  not  available  at  present  except  some  atterrpts  in  laminar  flows** '74,76. 
However,  some  calculations  will  appear  in  the  1980-01  Stanford  Conference^ . 

Although  the  origin  of  the  comer  flow77  j.s  Reynolds-shear-stress-driven  and 
is  different  from  the  jucntion  flow  considered  here,  the  basic  flow  features  have 
some  similarities  and  the  numerical  calculations78,79  of  the  comer  flow  can 
provide  helpful  information. 

If  momentum  integral  methods  are  attenpted,  then  the  mean  velocity  profiles, 
particularly  in  the  crosswise  direction,  must  be  constructed  based  on  experimental 
data  such  as  the  secondary  flow  vectors  shewn  in  Figure  lib.  No  such  attempt  has 
been  made  at  present.  According  to  the  experimental  evidenced,  the  strength  of  the 
streamwise  vortex  decreases  and  the  size  increases  downstream  due  to  the  diffusion  of 
the  Reynolds  shear  stresses.  This  diffusion  effect  and  the  mean  velocity  profiles 
are  the  main  features  that  need  to  be  considered  in  the  modeling  of  the  junction 
flow. 

Relatively  abundant  and  successful  calculations  have  been  performed  in  the 
flow  fields  relevant  to  region  (3)  where  the  non-crossover  crossflow  is  important. 
However,  when  crossover  of  the  crossflow  occurs,  numerical  calculations  encounter 
difficulties  and  in  general  the  results  are  not  satisf actory55 .  a  review  of  the 
numerical  methods  by  finite  differential  approaches  with  emphasis  on  turbomachinery 
has  been  given  by  McNally  and  SockolSO.  A  parallel  review  based  on  momentum 
integral  approaches  has  been  given  by  Horlock  and  Perkins**!. 
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Because  of  its  simplicity,  the  momentum  integral  approach  is  attractive  in 
this  flow  region.  The  most  popular  crossflow  models  are  those  by  Mager  and  Johnston 
and  those  using  polynomials.  For  a  general  discussion  of  these  models,  one  is 
referred  to  Reference  11.  A  comparison  of  the  results  obtained  by  using  the 
Mager  and  the  Johnston  model  has  been  made  by  Myring82.  it  is  found  that 
because  the  Johnston's  model  includes  a  Reynolds  number  dependency,  it  leads  to 
a  better  agreement  with  the  experiment  than  the  Mager' s  model  at  high  Reynolds 
numbers.  However,  comparisons  made  by  9nith83  on  several  experiments  have  not 
indicated  superiority  of  either  model.  Because  of  the  rather  artificial 
requirement  of  satisfying  the  boundary  conditions  at  the  wall  and  the  edge  of 
the  boundary  layer  in  order  to  determine  the  coefficients  of  the  polynominal, 
the  polynomial  model  does  not  appear  to  be  general.  At  present,  the  attenpts  to 
include  the  S-shaped  cross  flew  have  only  produced  a  qualitative  agreement  with 
the  experiment  at  best.  Some  discussions  of  this  subject  are  given  in 
References  11  and  84-86. 

WAKE 

It  is  known  that  Reynolds  shear  stresses  play  important  roles  in  the  rapidly 
relaxing  turbulent  flow  in  near  wakes.  A  review  of  two-dimensional  experimental 
data  and  prediction  methods  has  been  given  by  Patel  and  Scheuerer**7.  They  also 
evaluated  the  use  of  the  k- <  turbulence  closure  model,  and  found  that  while  the 
results  in  near  wakes  are  reasonable,  the  results  in  far  wakes  are  not  in  satis¬ 
factory  agreement  with  experimental  data.  The  techniques  used  in  the  prediction 
of  near  wakes  of  airfoils  are  a  good  source  to  search  for  a  suitable  method.  Some 
current  reviews  of  the  subject  are  given  by  Lock88,  Look  and  Firmin89,  Melnik90 
and  Le  Balleur  et  al.91 


Sene  simpler  approaches  to  predict  near  wakes  are  of  interest.  For  a  body 
of  revolution,  Wang  and  Huang 34  use  a  polynomial  to  connect  the  stem  boundary 
layer  region  and  the  far  wake  region  while  Nakayama  et  al58  enforce  the  conservation 
of  momentum  in  the  longitudinal  direction.  The  treatment  of  wakes  using  a  momentum 
integral  method  has  been  discussed  by  Green92,  and  subsequently  improved  by  Green 
et  al65.  in  this  method,  the  governing  equations  in  the  wake  are  the  same  as  in  the 
boundary  layer  except  that  the  skin  friction  is  set  equal  to  zero  since  there  is 
no  wall  boundary  and  the  dissipation  length  scale  is  doubled  in  order  to  match  the 
behavior  in  the  far  wake. 


SUMMARY  AND  RECOMMENDATION 

Due  to  the  limitation  of  computer  capability,  solving  the  time-dependent 
full  Navier-Stokes  equations  for  the  purpose  of  the  wake  prediction  of  submerged 
appended  bodies  is  not  feasible  at  present.  Even  when  time  dependency  is  eliminated 
by  the  time-averaging  method,  the  governing  equations  are  of  the  elliptic  type  and 
the  computing  effort  for  the  problem  considered  here  is  still  too  enormous  to  be 
practical.  Therefore,  further  simplifications  must  be  made.  For  example,  the 
diffusion  of  momentum  in  the  predominant  flow  direction  can  be  neglected  such  that 
the  elliptic  type  is  simplified  to  the  partially  parabolic  type;  or  the  thin 
boundary  layer  approximation  supplemented  by  viscous/inviscid  interaction  can  be 
made.  Once  the  type  of  the  governing  equations,  the  turbulence  closure  model  and 
the  calculation  method  of  the  potential  flow  have  been  chosen,  the  numerical  method 
can  proceed  according  to  either  the  finite  differential  approach  or  the  momentum 
integral  approach. 

If  the  finite  differential  approach  is  followed,  then  the  partially  parabolic 
type  of  the  governing  equations  and  a  turbulence  closure  model  no  more  carpi icated 


than  the  k-(  model  are  preferred.  Numerical  experiments  indicate  that  the 
partially  parabolic  type  of  the  governing  equations  supplemented  by  the  k-< 
model  will  likely  be  adequate  for  the  flow  regions  dominated  by  the  streanwise 
vortex  motion  (region(2))  and  by  the  crossflow  (region  (3)).  However,  its 
adequacy  in  the  region  where  flow  reversal  in  the  predominant  flow  direction 
occurs,  i.e.  region  (1)  where  the  horseshoe  vortex  is  generated,  is  uncertain. 
If  the  latter  method  is  adopted,  then  the  generality  of  the  finite  differential 
approach  will  be  limited  to  the  class  of  problems  where  the  experimental 
data  used  are  applicable. 

The  major  shortcoming  in  the  use  of  the  momentum  integral  approach  to 
threedimensional  flow  problems  is  the  fact  that  the  mean  velocity  profile 
in  the  crossflow  direction  varies  significantly  from  one  flow  problem  to 
the  other,  and  an  adequate  similarity  law  of  the  mean  velocity  profile  in 
the  crossflow  direction  applicable  to  a  large  class  of  flow  problems  has  yet 
to  be  found.  However,  if  the  aim  of  the  application  is  not  universal,  but 
is  confined  to  the  special  class  of  flows  around  appended  bodies  of  limited 
geometrical  variations,  then  the  more  empirically-oriented  momentum  integral 
approach  will  most  likely  be  adequate.  The  continued  usefulness  of  momentum 
integral  approach  to  three-dimensional  flow  problems  has  been  demonstrated 
in  seveal  recent  international  conferences  where  finite  differential  methods 
and  momentum  integral  methods  were  compared  and  evaluated.  If  the  momentum 
integral  approach  is  adopted  for  the  wake  prediction  of  submerged  appended 
bodies,  then  the  main  task  will  be  to  construct  the  empirically  based  models 
of  mean  velocity  profiles  appropriate  for  the  three  flew  regions 
mentioned  above. 
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The  choice  of  a  numerical  method  is  in  general  based  on  the  consideration  of 
generality  versus  practicality  and  accuracy  versus  efficiency.  There  is  a  common 
agreement  that  finite  differential  methods  tend  to  be  more  general  and  to  provide 
more  accurate  and  detailed  solutions,  or  have  the  potential  to  do  so;  while  momentum 
integral  methods  are  less  general,  and  are  only  capable  of  providing  solutions  of 
global  flow  features.  On  the  other  hand,  momentum  integral  methods  are  much  simpler; 
therefore  are  conputationally  more  efficient  and  are  practical  for  many  engineering 
applications.  It  has  been  mentioned  above  that  the  elliptic  type  of  governing 
equation  is  adequate  for  the  problem  considered  here,  but  the  computational  cost 
and  effort  are  prohibitive  at  present.  If  the  simpler  partially  parabolic  type  of 
governing  equation  is  adopted,  then  it  is  likely  that  the  computational  difficulties 
in  thick  turbulent  boundary  layers  and  large  cross  flows  can  be  overcome.  However, 
the  capability  of  the  partially  parabolic  type  equations  in  computing  the  flew 
region  where  the  horseshoe  vortex  is  generated  has  not  been  demonstrated. 

From  the  above  discussions,  it  appears  that  at  the  present  time,  there  is  no 
adequate  method  which  is  capable  of  computing  the  wakes  of  submerged  appended  bodies 
in  a  complete  fashion,  and  at  the  same  time,  meeting  the  criteria  of  both 
generality  and  accuracy  of  the  solution.  For  this  reason,  two  approaches  are 
recommended  to  be  pursued  simultaneously.  One  approach  is  for  immediate  applications 
and  the  other  approach  for  long-range  research.  For  the  first  approach,  the  momentum 
integral  method  is  recommended  due  to  its  relative  simplicity  and  its  flexibiilty 
to  adapt  to  experimental  data.  For  the  second  approach,  a  numerical  procedure  no 
less  sophisticated  than  the  one  by  solving  the  partially  p>arabolic  Navier-Stokes 
equations  with  the  K-f  turbulence  model  should  be  considered. 
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Figure  1  -  Variation  of  Boundary-Layer  Thickness  and  Local 
Radius  of  the  Body  (Taken  from  Reference  1) 


Figure  2  -  Variation  of  Wall  Static  Pressure  and  Total  Pressure 
at  the  Edge  of  the  Boundary  Layer. 
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(Taken  from  Reference  1) 


Figure  5  -  Flow  Field  of  the  Appendage/Body 
Junction  (Taken  from  Reference  5) 


Figure  6  -  Wake  of  the  Appendage/Body  Junction  Induced 
Flow  Field  (Taken  from  Reference  5) 
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Figure  7  -  Non-Crossover  Type  of  Crossflow 
(Taken  from  Reference  11) 


Figure  8  -  Crossover  (or  S-shaped)  Type  of  Crossflow 
(Taken  from  Reference  11) 
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SECONOARY  SEPARATION  LINE 
LAMINAR  / 


Figure  10  -  Comparison  Between  Experimental  Laminar  and  Turbulent  Flows  About 
2:1  and  4:1  Ovals  (Taken  from  Reference  10) 


figure  11a  —  Idealized  Wing/Body  Junction 
(Taken  from  Reference  13) 


Figure  lib  -  Secondary  Flow  Vectors  at  a  Downstream  Station 
X=1223  mm  of  the  Idealized  Wing/Body  Junction 
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EXPERIMENT 


TABLE  1 


Levels  of  Approximation  and  the  State-of-the-Art  of 
Computational  Fluid  Dynamics 
(Taken  from  Reference  22,  1980) 


Level  of  Approximation 


State-of-the-Art 


1  -  Ideal  fluid 

2a  -  Ideal  fluid  +  nondimensional  enpirical  correlations 


2b  -  Ideal  fluid  +  boundary  layers  (weak  coupling)  x 

2c  -  Parabolized  approximations  of  the  Navier-Stokes  x 

equations 

2d  -  Ideal  fluid  +  isobaric  separated  zones  x 

3a  -  Ideal  fluid  +  thin  viscous  layers  (strong  coupling)  x 

3b  -  Ideal  fluid  +  Navier-Stokes  zones  (strong  coupling) 


4  -  Direct  overall  solutions  of  the  Navier-Stokes  equations: 

4a  -  Averaged  equations  +  turbulence  model 

4b  -  Filtered  equations  (simulation  of  the  large  turbulent 
structures ) 

4c  -  Complete  equations  (numerical  simulation  of  turbulence). 
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Development  of  Computational  Aerodynamics 
(Taken  from  Reference  24,  1979) 


TABLE  4 


Type  of  Equations 

Hyperbolic 30 

Parabolic31 

Partially 

parabolic32 

Elliptic33 


Classification  of  the  Governing 


Approximations 

Thin  boundary  layer 
and  no  diffusion 

Thin  boundary  layer 

Diffusion  in  the  stream- 
wise  direction  neglected 

No  or  few  approximatiofis 


Equations 


Flow  Problems  Solved 


Flews  with  a  predominant 
direction;  no  vortex 
motion 


Same  as  above 


Curved  duct  flows. 
Seme  degree  of  thick 
boundary  layer 


Horseshoe  vortex  and 
thick  boundary  layer 


TABLE  5 


Comparison  of  Nunerical  Methods 


Major 

Editor 

(year) 

Type  of  Problem 

Numerical  Methods 

Differen¬ 
tial  Integral 

Comments 

Kline38 

(1968) 

33  test  cases; 
mostly  2-D 
turbulent 
boundary  layers 

9 

20 

Of  7  in  the  top  class, 

4  are  momentum  integral 
methods  and  3  are  finite 
differential  methods. 

East52 

(1975) 

7  test  cases, 
some  with 
separation 

6 

3 

Momentum  Integral  methods 
have  impressive,  if  not 
better  performances. 

Humphrey^ 

(1978) 

A  3-D  wing 

5 

3 

Despite  the  fact  that  the 
flow  is  not  strongly  3-D, 
the  computed  results 
scatter  by  as  much  as  20%; 
no  preference  to  either 
numerical  method  indicated. 

Lindhout54 

(1979) 

A  3-D  wing 
with  separation 
in  the  root 
section 

5 

4 

The  potential  capability 
of  finite  differential 
methods  has  not  been  demon¬ 
strated  .  Momentum  integral 
methods  appear  to  provide 
more  consistent  results. 

Larson55 

(1980) 

2  test  cases  of 

3-D  Ship 
boundary  layers 

7 

10 

Both  finite  differential 
and  momentum  integral 
methods  failed  to  obtain 
adequate  solutions  in  the 
thick  turbulent  boundary 
layer  near  the  stem. 

Kline56 

(1980-81) 

66  test  cases  of 
computer 
turbulent  flows 

A  total  of  73 
different  methods 

Both  finite  differential 
and  momentum  integral 
methods  are  well  worthy 
of  further  study  and 
refinement. 
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TABLE  6* 


Methods  to 


Methods 

1.  Elliptic 
equations 

2.  Partially 
parabolic 
equations 


3 .  Presure 
iteration 

4 .  Viscous/inviscid 
matching  at  the 
height  of  the 
boundary  layer 


5.  Viscous/inviscid 
matching  at  the 
height  of  the 
displacement 
thickness 

6 .  Vortex  Method 


7 .  Strong  surface 
transpiration 
method 


*The  order  of  listing 
computing  efforts. 


Account  for  Thick  Turbulent  Boundary  Layers 


Comments 


References 


Computer  cost  prohibitive;  Phoenics  Code33 

not  practical 

Iterative  marching  integra-  Abdelmequid  et  al35 

tion  scheme.  The  rate  of 

convergence  of  solution 

depends  strongy  on  the 

quality  of  the  initial  guess 

on  pressure 

2-D  method  has  been  developed  Mahgoub  &  Bradshaw^? 


An  application  to  axisynmetric  Nakayama  et  al58 

bodies  of  revolution  has  been 

considered 


Less  accurate  than  -he  above  Wange  &  Huang59 

method 


The  tail  region  of  the  body  of  Geller^O 
revolution  is  calculated  by 
the  vorticity  equation 

The  pressure  of  the  real  flow  is  Le  Balleur^l 

set  equal  to  that  of  the  in- 

viscid  flow  inside  the  boundary 

layer  in  the  first-order 

approximation 


is  according  to  the  decreasing  degree  of  demand  on 
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fit  the  data;  Johnston's  model  has  a 
better  fit 

56. 

X 

Turbulent  flow  about  a 
teardrop  body  mounted  on 
a  flat  plate 

Only  the  first 
quadrant  of  the 
leading  edge  region 
is  considered 

The  flow  in  the  separated  region 
is  character i zed  by  a  single  con¬ 
centration  of  vorticity  (horseshoe 
vortex) 

57. 

X 

X 

2-D  separate!  turbulent 
boundary  layers 

Lag-entrainment 

method 

An  exploratory  wrk  on  the  use  of 
momentum  integral  methods  to 
calculate  separated  flows 

58. 

X 

X 

Structure  of 
equilibrium  turbulent 
boundary  layers 

2-D 

Seven  turbulent  boundary  layers  rang¬ 
ing  from  mildly  favorable  to  severe 
adverse  pressure  gradients  were  shown 
to  be  good  approximations  to  2-D 
equilibria  flows 

59. 

X 

X 

X 

3-D  turbulent  boundary 
layer  approaching 
separation 

Thin  boundary  layer 

Bradshaw's  one-equation  turbulence 
model  is  used  but  it  fails  to  predict 
separation 

60. 

X 

X 

I 

X 

Review  of  experiments  and 
predictions  of  3-D 
turbulent  boundary  layers 

•  • 

No  computational  method  of  the  flow  in 
the  wing -body  junction  is  presented  in 
the  colloquium 

61. 

X 

X 

3-D  boundary  layer  and 
wake  over  a  wing/body 
ccrrbl  nation 

lag-entrainment 

equation 

The  horseshoe  vortex  generated  in  the 
wing/body  junction  is  suppressed 

62. 

X 

• 

Turbulent  horseshoe 
vortex  flow 

Thin  boundary  layer 

The  inverse  mode  of  computation  is 
used  to  march  across  the  separation 
line 

63. 

X 

a 

Thick  turbulent  boundary 
layer 

Body  of  revolution 

The  tail  region  is  calculated  by  the 
vorticity  equation 
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64. 

X 

X 

Turbulent  comer  flow 

Based  on  energy  and 

vorticity 

consideration 

The  origin  of  the  secondary  flow  in 
the  comers  is  due  to  the  transverse 
gradients  of  the  Reynolds  shear 
stress 

65. 

X 

X 

Turbulent  comer  flow 

Thin  boundary  layer 
approximation  with 
an  algebraic 

Reynolds  stress 
model 

Predicted  Reynolds  stress  distributions 
were  compared  with  experimental  data 

66. 

x 

X 

Turbulent  boundary 
layer  and  wake 

2-D 

The  main  features  are  the  treatment 
of  the  wake  and  the  formulas  for  skin 
friction  and  shape  factor  are 
derived  from  the  zero  pressure 
gradient  case 

67. 

X 

X 

Turbulent  boundary 
layer  and  wake 

2-D 

Derives  the  lag-entrainment  equation 
which  is  the  integral-method  version 
of  Bradshaw's  one-equation  turbulence 
model 

68. 

X 

X 

Secondary  flow  about 
struts  and  airfoils 

Inviscid  theory 

Comparisons  between  theory  and  experi¬ 
ment  are  made 

69. 

X 

X 

Turbulent  boundary 
layer 

2-D 

The  use  of  die  entrainment  equation 
as  an  auxiliary  equation  is  introduced 
for  the  first  time 

70. 

X 

Review  of  subgrid 
scale  modeling 

Prediction  methods  based  on  subgrid 
scale  modeling  are  important  for  flows 
where  the  range  of  length  scales  is 
very  wide 

71. 

X 

Potential  flow 

Panel  method 

Enphasis  is  on  the  quadra lateral  panel 
method 

72. 

X 

X 

Viscous  flows  over 
ship  stems 

Hun  boundary 
layer 

Hie  viscous-inviscid  interaction 
approach  is  strongly  criticized  due  to 
a  non-local  reaction  of  the  boundary 
layer 

73. 

X 

X 

X 

Thick  turbulent 
boundary  layers 

Body  of  revolution 

A  modified  displacement  thickness 
method  is  used  to  account  for  the 
normal  pressure  gradient 

74. 

X 

* 

Crossflow  models  in 
turbomachine 

Available  crossflow  models  are  reviewed 
and  a  new  model  is  proposed 

75. 

X 

X 

X 

X 

A  review  of  3-D 
turbulent  boundary 
layers 

The  prediction  of  separation,  comer 
flows  and  secondary  flows  are  discussed 

_  -  -  — 
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76. 

X 

■ 

Flows  in  turbcmachines 

Inviscid  theory 

The  concept  of  averaging  the  flow 
across  the  pitch  (i.e.  the  passage- 
averaged  flow)  is  introduced 

77. 

X 

X 

Review  of  secondary  flows 

Both  theory  and  experiment  are 
discussed  with  an  emphasis  on 
turbomachine  applications 

70. 

X 

B 

Review  of  flows  in 
turbcmchines 

The  application  of  momentum  integral 
methods  is  emphasized 

79. 

X 

Turbulent  flow  about  a 
cylinder  mounted  on  a 
flat  plate 

Only  the  first 
quadrant  in  the 
leading  edge  of 
the  cylinder  is 
considered 

Johns ton's  mean  velocity  profile  model 
lias  a  better  correlation  than  that  of 
Coles'  model 

00. 

X 

Thick  turbulent 
bourdary  layers 

Body  of  revolution 

Pressure  distributions,  mean  velocity 
profiles  and  Reynolds  shear  stresses  . 
are  measured.  A  parallel  work  to 

Bibl.  146  and  143 

81. 

X 

5 

3 

Comnarison  of  numerical 
methods  or.  a  3-D  wing 
turbulent  boundary 
layer 

1  test  case 

Despite  the  relative  sinplicity  ot  t'.e 
test  case,  the  over-all  prediction  is 
not  very  good 

G2. 

X 

X 

3-D  separation  of  flews 

Application  of  topology  to  flow 
visualization  of  the  flows  around 
free  or  surface-mounted  obstacles 

03. 

X 

Comparison  of  numerical 
methods 

3  test  cases 

The  flew  around  a  cylinder  mounted 
on  a  flat  plate  is  one  of  the  test 
cases 

84. 

X 

2-D  turbulent  boundary 
layer 

Coles'  model 

Instead  of  the  entrainment  equation, 
the  energy  equation  is  used  and  is  sh 
to  he  adequate  even  for  non-equilibrium 
turbulent  boundary  layers 

85. 

X 

Secondary  flow  of  3-D 
turbulent  boundary 
layers 

Hie  Johnston's  triangular  (mean 
velocity  profile)  model  was  bom 

86. 

X 

3-D  turbulent  boundary 
layer  with  a  plane 
of  symmetry 

Johnston's 

triangular 

model 

A  good  agreement  between  calculation 
and  measurement  is  obtained  in  the 
non-separated  region 

87. 

X 

l 

3-0  turbulent  boundary 
layer  induced  by  a  step 

_ 

Ttie  angle  between  the  mean  velocity 
gradient  vector  and  the  shear  stress 
vector  is  significant  even  close  to  the 
wall.  Thus,  the  eddy-viscosity  model 
is  not  adequate  for  the  flow  considered 
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08. 

X 

■ 

Turbulent  boundary 
layers 

2-D 

A  general  family  of  the  velocity 
profiles  of  2-D  movirg-equilibriLm 
turbulent  boundary  layer  was  derived 

09. 

X 

1 

I 

Discrete  vortex  metiiod 

2-D 

The  method  is  used  to  calculate  the 
unsteady ,  separated  flow  over  a  2-D 
airfoil 

90. 

X 

X 

1 

I 

Decay  of  streamwise 
vorticity 

The  streamwise  vorticity  is  found  to 
persist  more  than  100  protuberance 
diameters  downstream  both  in  laminar 
and  turbulent  bounday  layers 

91. 

X 

1 

1 

Comparison  of  nunerical 
metiiods  on  complex 
turbulent  flows 

66  test  cases 

The  emphasis  was  on  2-D  and  3-D 
complex  turbulent  flews.  Integral 
methods  continue  to  perform 
adequately. 

92. 

X 

9 

20 

Comparison  of  numerical 
methods  on  mainly  2-D 
turbulent  boundary 
layers 

33  test  cases 

No  distinct  superiority  or  finite 
differential  methods  over  integral 
methods  were  observed 

93. 

X 

UirtKilent  shear  flows 

The  current  status  of  prediction 
methods  and  ten-year  outlook  are 
discussed  (in  1978) 

94. 

X 

Turbulent  ixxindary  layer 
flow  in  a  S-shaped 
channel 

Internal  flows 

No  mean  velocity  profile  models  are 
adequate  for  the  resulting  crossover 
type  of  crossflow 

95.  | 

X 

Modeling  of  crossflows 

Models  of  mean  velocity  profiles 
include  crossover  crossflows  are  p 
proposed 

96. 

X 

Ship  boundary  layer 

.. 

A  general  discussion  on  the  effect  of 
surface  curvature  upon  the  calculation 
of  turbulent  boundary  layers  of  ship 
hulls 

97. 

X 

Secondary  flows 

The  governing  equations  of  secondary 
flews  are  derived  in  a  rotating 
coordinate  system 

98. 

X 

Rview  of  the 
characteristics  of 
ship  boundary 
layers 

Surface  ships 

The  use  of  integral  methods  to  ship 
boundary  layers  is  criticized 

99. 

X 

1 

X 

Review  of  numerical 
methods  on  ship 

Ixxindary  layers 

Surface  ships 

The  use  of  integral  methods  to  ship 
boundary  layers  is  criticized 

99. 


Review  of  numerical 
methods  on  ship 
Ixxindary  layers 


Surface  ships 


The  use  of  integral  methods  to  ship 
boundary  layers  is  criticized 
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Modeling  of  the  mean 
velocity  profile 


Turbulent  flow 
within  a  turbine 
cascade  passage 

10  Comparison  of  numerical 
methods  on  3-D  ship 
boundary  layer 


Review  of  turbulence 
closure  models 


Turbulent  flows 


Turbulent  boundary 
layers  on  spinning 
and  curved  surfaces 

Survey  of  numerical 
methods  on  viscous 
flows  at  high 
Reynolds  numbers 

Turbulent  boundary 
layers  and  wakes 
including  separation 


Viscous  flows 


Boundary  layer  on  3-D 
wing 


APPROXIMATIONS 
OR  MODELS 


Based  on  the  data 
obtained  by 
Lang ban  et  al 
(Bibl.  101) 


2  test  cases 


High  Reynolds 
numbers 


No  crossflow 


Review  of  vortex  methods  I  2-D  and  3-p# 


A  model  of  crossflow  including  the 
crossover  type  is  proposed 


In  addition  to  flow  visualization, 
both  static  pressure  and  mean 
velocity  are  measured 

Both  integral  and  differential 
methods  have  difficulties  in  the 
stem  region  where  the  turbulent 
boundary  layer  is  thick 

The  emphasis  is  on  two-equation  or 
stress-equation  models 

The  reasons  for  choosing  the  K-E 
(two-elation)  turbulence  model  are 
carefully  explained 

The  effects  of  surface  curvatures  are 
investigated  in  view  of  modifying  the 
turbulence  model 

The  state-of-the-art  is  examined 


A  new  concept  of  visoous/inviscid 
interation  and  matching  is  proposed. 
In  this  method ,  both  viscous  and 
inviscid  flows  are  calculated 
simultaneously 

Nimerical  methods  of  both  finite 
differential  and  momentum  integral 
methods  are  discussed 

The  surface  transpiration  method  is 
found  to  be  equivalent  to  the 
displacement  surface  method  in  adding 
viscous  corrections 

Vortex  methods  are  potentially 
powerful  methods  to  predict  high- 
Reynolds  nunber  flows  particularly 
when  separation  occurs 
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X 

Turbulent  boundary 
layers 

2-D 

Comparison  of  the  use  of  the  entrain¬ 
ment  energy  and  moment  of  momentum 
equation  as  an  axiliary  equation 

Displacement  thickness 

2-D  and  3-D 

The  equivalence  of  the  displacement 
surface  method  and  the  surface  tran¬ 
spiration  method  is  illustrated 

Separation  in  flows 

A  good  review  of  the  fundamental 
structure  and  definition  of  3-D 
separation 

4 

Comparison  of  numerical 
methods 

Root  section  of  a 
3-D  wing 

Momentum  integral  metlxxis  gave  more 
consistent  results 

3-D  boundary  layers  on 
wings 

Simple  eddy- 
visoosity 
turbulence  model 

Despite  the  anisotropic  natut !  of  the 
turbulent  flow,  the  simple  edly- 
viscosity  model  gives  season  x>le 
results 

Review  of  niroerical 
methods  on  boundary 
layers  and  wakes  of 
airfoils 

Mostly  2-D 

Integral  methods  are  f-pnasized 

Review  of  numerical 
methods  used  in 
external  aerodynamics 

2-D  and  3-D 

New  developments  in  integral  metlxxis 
are  emphasized 

Viscous  secondary  flow 

Secondary  flow  is 
considered  as  a 
small  perturbation 
to  a  3-D  flow 

The  inclusion  of  the  diffusion  of 
vorticity  improves  the  agreement  with 
the  data 

Review  of  the  numerical 
modeling  of  vortex  flows 

The  appliation  to  ship  hydrodynamics 
is  emphasized 

Review  of  the  numerical 
modeling  of  vortex  flows 

A  symmetric  "winy" 
mounted  on  a  flat 
plate 

TVo  mean  velocity  components  and  three 
Reynolds  shear  stress  components  are 
measured  in  the  wing/body  junction 

Review  of  computation 
methods  used  in 
turtonachine 

A  thorough  review  of  various  finite 
differential  methods 

• 

Thick  boundary  layer 

2-D 

The  normal  pressure  gradient  is 
calculated  by  an  iterative  schej.ie 
without  actually  solving  the  elliptic 
system  of  equations 
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123. 

X 

X 

3-D  laminar  boundary 
layer 

Moment  of 
monenbm  equation 

A  method  for  stable  transition  frcm 

2-D  computation  to  3-0  confutation 
is  discussed 

124. 

X 

X 

3-D  turbulent  boundary 
layers  on  ship  hulls 

two-equation 
model  turbulence 

A  further  exposition  of  Spalding's 
partially  parabolic  calculation 
procedure 

125. 

X 

X 

Flow  around  an  appended 
submerged  body 

two-equation 
turbulence  model 

A  sinplified  version  of  the  3-D 
partially  parabolic  scheme  is  used. 

Ttie  prediction  of  the  horseshoe 
vortex  is  not  attepted 

126. 

X 

X 

Flows  in  turbcmachines 

passage-averaged 

flow 

The  governing  equations  for  wail 
boundary  layers  in  turbcmachines  are 
derived 

127. 

X 

Equilibrium 
turbulent  boundary 
layers 

2-D 

The  emphasis  is  on  the  velocity  defect 
profiles  (i.e.  outer  region  of  the 
boundary  layer) 

128. 

X 

X 

End-wall  boundary 
layer  theory  in 
turbomachines 

3-D 

Seme  approximate  profiles  of  the 
streamwise  and  the  crossflow  velocity 
are  assumed 

129. 

X 

Review  o f  turbulence 
closure  models 

Review  includes  one-equation,  two- 
equation  and  stress-equation  models 

130. 

X 

Review  of  numerical 
metliods  on  boundary 
layers  and  wakes  of 
airfoils 

2-D 

Both  weak  and  strong  viscid-inviscid 
interactions  are  discussed.  The 
asymptotic  expansion  method  is 
emphasized 

131. 

X 

Turbulent  boundary  layer 
over  a  longitudinally 
curved  surface 

A  data  base  for  the  flow  over  concave 
and  convex  surfaces  is  established. 

The  comparison  with  calculations 
based  on  Bradshaw's  technique  (Bibl. 

8  and  9)  gives  good  agreement 

132. 

X 

X 

Thick  boundary  layer 
and  wake 

Partially  parabolii 
and  two-equation 
turbulence  model 

The  superiority  of  the  two-equation 
turbulence  model  is  claimed 

133. 

X 

X 

3-0  turbulent  boundary 
layer 

Non-orthogonal 
coordinate  system 
and  entrainment 
equation 

Discussions  of  the  characteristics 
and  the  miner ical  procedures  of  the 
governing  equations  are  given 

134. 

X 

X 

Boundary  layer  and  wake 
of  a  body  of  revolution 

Coles'  model  and 

entrainment 

equation 

The  displacement  surface  method  is 
used  to  account  for  the  normal 
press  gradient  in  the  thick  turbulent 
boundary  layer 
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■ 

X 

X 

3-D  turbulent  boundafy 
layer  approaching 
separation 

Hager's  crossflow 
model  and 
entrainment 
equation 

Interest  centers  on  the  prediction  of 
the  pressure  rise  to  separation 

136. 

X 

X 

Ttii.dk  turbulent  boundary 
layer  and  near  wake 

A  linear  crossflow 
profile  assumed 

Additional  terms  due  to  surface 
curvature  effects  are  included 

137. 

X 

X 

X 

Review  of  3-D  turbulent 
ooundary  layer 

Detailed  derivations  of  the  governing 
equations  are  given.  Qnphasis  are  on 
momentum  integral  metliods 

138. 

X 

X 

X 

Review  of  advances  in 

3-D  turbulent  boundary 
layer 

Both  finite  differential  and  momentum 
integral  methods  are  surveyed 

139. 

X 

Missiles  aerodynamics 

Inviscid  rotational 
flow 

The  forces  on  wing/body  combinations 
at  angles  of  attack 

140. 

X 

X 

X 

Ship-stern  bounday  layer 

Surface  ships 

Contains  an  extensive  literature  survey 

141. 

X 

Turbulent  flow  in  a  wing/ 
body  junction 

A  symmetric  "wing" 
mounted  on  a  flat 
plate 

Pressure,  mean  velocity  coiponents  and 
RMS  velocity  fluctuation  are  measured 
at  two  boundary  layer  thicknesses 

142. 

X 

Review  of  spectral  methods 

The  use  of  spectral  methods  for  the 
prediction  of  viscous  flow  is  still  in 
the  research  stage 

143. 

X 

X 

X 

Review  of  flews  in 
turbomachines 

3-D 

Various  experimental  data,  theories 
and  computation  methods  of  end-wall 
boundary  layers  in  turbomachines 
are  discussed 

144. 

X 

X 

Calculation  procedure  for 
parabolic  flows 

3-D 

The  pressure  in  the  streamwise 
momentum  equation  is  decoupled  by 
averaging  over  the  cross-section. 

Then  a  calculation  procedure  with 
forward  marching  method  is 
established. 

145. 

X 

X 

X 

Thiele  turbulent  boundary 
layer 

Body  of  revolution 

Equations  adequate  for  thick 
turbulent  boundary  layers  are  derived, 
but  the  effect  of  the  longitudinal 
curvature  is  not  included. 

146. 

X 

• 

Thick  turbulent  boundary 
layer 

Modified  spheroid 

Pressure  distributions,  mean  velocity 
profiles  and  Reynolds  shear  stresses 
are  measured 

147. 

X 

Thick  turbulent  boundary 
layer  and  near  wake 

bsw-drag  body  of 
revolution 

Pressure  distributions,  mean  velocity 
profiles  and  Reynolds  shear  stresses 
are  measured 
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Wear  and  Car  wake 

2-D 

The  model  is  shown  to  give  good 
results  in  near  wakes,  but,  does  not 
give  correct  results  in  far  wakes. 

Thick  turbulent  boundary 
layer  and  wake 

Body  of  revolution 

The  calculations  by  differential  and 
integral  methods  are  compared  with  the 
measurements.  No  significant 
difference  in  accuracy  between  the  two 
methods  is  found. 

3-D  laminar  and  turbulent 
boundary  layers 

Body  of  revolution 
at  incidence 

The  need  for  a  more  a<3equate  strong 
visoous/inviscid  interaction  scheme  is 
indicated  for  this  type  of  flow. 

Laminar  Clows  about 
circular  cylinders 
mounted  on  a  flat  plate 

Small  crossflow 

Agreement  between  calculation  and 
measurement  is  in  general  good  before 
separation.  The  inaccuracy  is  blamed 
on  inadequate  potential  flow 
calculation. 

Rankine  ovals  mounted  on 
a  flat  plate 

Small  cross  flow 

Same  raiments  as  in  Bibi.  151. 

The  separated  regions  are  smaller  when 
the  boundary  layers  are  turbulent. 

Review  of  3-D  inter¬ 
actions  and  vortical 
flows 

Mostly  inviscid 

The  emphasis  is  on  high-speed 
aircraft.  An  extensive  list  oE 
references  is  provided. 

Turbulent  comer  flow 

The  mean  streanwise  vorticity  is  shown 
to  arise  from  mean  flow  skewing  and 
the  inhemogeneity  of  anisotropic  wall 
turbulence. 

Modeling  of  the  mean 
velocity  profile 

Based  on  the  data 
of  Homung  & 

Jonbert 

The  inner  and  outer  regions  of  the 

3-D  turbulent  boundary  layer  ace 
considered. 

General  purpose  viscous 
flew  aodes 

A  computer  code  for  solving 
parabolic,  hyperbolic  and  elliptic 
type  of  tiie  governing  equations. 

Turbulent  flow  in  the 
leading  edge  region  of  a 
teardrop  body  on  a  flat 
plate 

Only  the  first 
quadrant  of ( the 
leading  edge  region 
outside  of  the 
separation  line  are 
considered 

Mean  velocity  component,  static 
pressure  and  skin  friction  are 
measured;  no  adequate  3-D 
similarity  law  is  found 
to  fit  the  data. 

Mean  velocity  profile  of 
the  flow  about  a 
teardrop  body  mounted 
on  a  flat  plate 

Data  further  downstream  of  those 
obtained  in  Bibl.  56  and  79  are 
measured. 
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K  8 


APPROXIMATIONS 
OR  MODELS 


Flow  In  curved  ducts 


Review  of  turbulence 
closure  models 


Review  of  both  turbulence 
closure  models  and 
numerical  methods 

Thick  boundary  layers 


Turbulence  closure  models 
for  3-D  boundary  layers 

Structure  of  turbulent 
bound ry  layers 

Separated  flows  due  to 
protuberances 


Wakes  of  tail-fin/body 
intersections 


Turbulent  flow  in  a  wing/ 
body  junction 


Modeling  of  crossflow 


Crossflows  in  3-0 
turbulent  bounday 
layers 


3-D  turbulent  boundary 
layer 


Turbulent  comer 
flow 


Flow  visualization 


Body  of  revolution 


A  symmetric  "wing" 
mounted  on  a  flat 
plate 

Based  on  a  poly¬ 
nomial  representa¬ 
tion 


The  partially  parabolic  system  with 
the  K-E  (two-equation)  turbulence 
model  is  applied  to  compute  duct  flows 

The  review  includes  models  of  zero- 
equation,  one-equation,  two-equation 
stress-equation  and  large  eddy 
simulation 

The  current  status  of  mcmentun 
integral  methods  is  also  discussed 


The  second-order  boundary  layer 
equations  are  derived  to  account  for 
effects  due  to  surface  curvatures 

Attempts  are  made  to  account  for 
anisotropic  shear  stress 

An  extension  of  the  previous  analysis 
Bibl.  156 

Multiple  vortices  are  observed.  A 
some  what  weak  correlation  between 
the  separation  distance  and  the 
protuberance  is  also  observed 

Demonstrates  the  importance  of  the 
interference  effects  of  tail-fin/ 
body  intersections 

Static  pressure,  skir.  friction,  mean 
velocity  components  and  Reynolds  shear 
streres  are  measured  along  the  junction 

The  (pi,qi)  family  of  crossflow  models 
are  proposed.  The  models  include 
cross-over  type  crossflow 


Thin  boundary  layer  Comparison  is  made  with  the  data  of 
approxnvation  with  skewed  turbulent  boundary  layers 
the  (pi,qi)  family 
of  crossflow  models 

Entrainment  Extension  of  Head's  entrainment 

equation  equation  for  incompressible  flow 

to  compressible  flow 

3-D  A  two  layer  eddy-viscosity  turbulence 

model  is  used 


TYPE  OF  PROBLEMS 

APPROXIMATIONS 

OR  MODELS 

COMMENTS 

Review  of  finite-element 
methods  in  fluid 
mechanics 

Only  fundamental  aspects  are 
discussed 

Review  of  separation 
phenomena  in  turbulent 
flow 

The  experimental  aspect  of  separation 
was  emphasized,  but  some  descriptions 
of  modelings  and  prediction  methods 
are  also  discussed 

Structure  of  separating 

2-D  turbulent  boundary 
layer 

Airfoil-type 

pressure 

distribution 

The  main  features  of  the  flow  in  the 
vicinity  of  separation  are  discussed 

Review  of  separation  in 
steady,  flows 

3-D 

Bnphasis  is  on  modeling  and  structure 
of  separation 

3-D  turbulent  boundary 
layer  on  an  infinite 
yawed  wing 

Entrainment 

equation 

The  results  obtained  by  using 
different  skin  friction  formulas 
and  mean  velocity  profiles  etc  are 
compared.  The  polynomial  type  of 
velocity  profiles  is  found  to  be 
unsatisfactory 

3-0  turbulent  boundary 
layers 

Non-orthogonal 
coordinate 
system  and 
entrainment 
equation 

Johnston's  crossflow  profile  is  found 
more  adequate  than  Mager's  profile  in 
the  problems  where  the  Reynolds  number 
is  high 

Secondary  flow  in  a 
cascade  passage 

Inviscid  theory 

Pioneer  work  in  secondary  flow 

Rectangular-wing /body 
and  delta-wing/body 
combinations 

Oil-flow  visualiza¬ 
tion  and  surface 
pressure  distribu¬ 
tions 

The  effects  of  flow  incidence  and  the 
Reynolds  nunber  on  the  strength  and 
the  position  of  the  vortices  ^re 
also  investigated 

Review  of  the  structure 
of  laminar  boundary 
layers 

The  triple-deck  structure  of  laminar 
boundary  layer  is  described 

Prediction  of  separation 
lines  and  calcuatlon  of 
separated  flows 

Inclined  ellipsoids 
and  infinite  swept- 
wlngs 

A  new  crossflow  profile  derived  from 
the  profile  of  a  flat  plate  is  used 

Prediction  of  separation 
in  turbulent  boundary 
layers 

2-0 

One  of  the  best  engineering  prediction 
methods  for  separation  in  2-D 
turbulent  boundary  layers 

Turbulent  boundary 
layers 

2-D 

A  new  family  of  mean  velocity  profiles 
is  proposed.  The  skin  friction 
obtained  from  this  family  is  an 
improvement  over  that  by  Ludwieg  and 
Tillmann 

Turbulent  boundary 
layers 


n 

NUMERICAL 

METHODS 

BIBLIOC.  NO. 

5 

i 

EXPERIMENT 

DIFFERENTIAL 

INTEGRAL 

TYPE  OF  PROBLEMS 

APPROXIMATIONS 

OR  MODELS 

COMMENTS 

184. 

X 

X 

Separations  in  Clows 

2-D  and  3-D 

The  topology  of  separations  occurring 
in  aircrafts  is  discussed 

X85. 

X 

Turbulent  boundary  layer 
near  separation 

The  behavior  of  the  turbulent  boundary 
layer  is  discussed;  a  fundamental  paper 

186. 

X 

Equilibria^  turbulent 
boundary  layers 

The  type  and  the  structure  of 
equal ibr ion  turbulent  boundary  layers 
are  discussed;  a  fundamental  paper 

187. 

X 

Law  of  the  wall 

3-D 

The  3-D  law  of  the  wall  was  derived 
with  the  effects  of  both  pressure 
gradients  and  inertial  forces  accounted 
for 

188. 

X 

X 

X 

3-D  turbulent  boundary 
layer  over  an  infinite 
swept -wing 

No  surface 

curvature 

effect 

Despite  the  anisotropy  of  the  shear 
stress,  computation  based  on 

Bradshaw's  model  give  good  results 

189. 

X 

X 

Separated  flows 

Laminar  flows 

The  quasi-simultaneous  method  is 
developed  as  an  alternative  of  the 
direct  and  the  inverse  mode 

190. 

X 

Separated  flows 

laminar  flows 

The  (juasi -simultaneous  method  is 
developed  as  an  alternative  of  the 
direct  and  the  inverse  mode 

191. 

X 

X 

Turbulent  boundry  layer 
in  a  ship  hull 

Thin  boundary 
layer 

A  new  crossflow  model  is  used 

192. 

X 

Classification  of  3-D 
boundary  layer  equations 

The  zones  of  influence  and  dependence 
are  discussed  in  terms  of 
subcharacteristics 

193. 

X 

X 

Flow  around  axisymnetric 
bodies 

Thin  boundary 
layer 

The  norrral  pressure  gradient  is 
accounted  for  by  the  displacement 
thickness  method 

194. 

X 

Potential  flow 

Panel  method 

The  triangular  panel  method  with  a 
linear  source  strength  is  developed 

195. 

X 

X 

Assessment  of  3-D 
turbulent  boundary 
layer  predicitcn 
methods 

10  test  cases 

The  test  cases  include  East  and  Hoxey's 
wing-body  junction  flow.  Bradshaw's 
rate-equation  (or  one-equation) 
turbulent*  model  is  recontnended 

196. 

X 

• 

Model  mg  of  mean 
velocity  profile 

2~D  flow 

An  analytic  expression  of  the  mean 
velocity  is  proposed  and  a  comparison 
with  the  measurement  is  also  made 

197 

X 

X 

Turbulent  boundary 
layers  with  separation 

Mean- flow  kinetic 
energy  integral 
{ instead  of  en¬ 
trainment 
equation) 

A  new  multi-layer  velocity  profile  is 
used 

84 
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